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Abstract: We study the Central Exclusive Production (CEP) of Xcj and rj c mesons 
at RHIC in proton-proton collisions. We consider the Xcj — >• J/ip + 7 decay channels 
and, recalling that the J z = suppression of the J = 1,2 states can be compensated 
by their larger branchings to J/tp + 7, present predictions of rates and distributions for 
Xc(o,i,2) production. Particular attention is paid to the impact of p± cuts applied to the 
outgoing protons, which can influence the relative rates significantly The distribution in 
the azimuthal angle difference between the outgoing protons and the proton p\ is also 
studied, and shown to depend sensitively on the spin and parity of the centrally produced 
meson, as well as being affected by the soft survival factors, S 2 . Two- and four-body decays, 
which are particularly relevant for XcO production, are also considered. We show that in 
the two-body case, backgrounds from 'direct' QCD production based on both perturbative 
and non-perturbative models are expected to be under control. 
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1. Introduction 

Central exclusive production (CEP) processes in high-energy collisions represent a very 
promising and novel way to study both QCD and new physics at hadron-hadron colliders. 
In particular, these reactions provide a valuable tool to investigate in detail the properties 
of resonance states, from 'old' SM mesons to BSM Higgs bosons (see for instance [l]-[6]). 
The CEP of an object X may be written in the form 

pp{p) ->• P + X + p{p) 

where '+' signs are used to denote the presence of large rapidity gaps. An attractive 
advantage of these reactions is that they provide an especially clean environment in which 
to measure the nature and quantum numbers (in particular, the spin and parity) of the 
centrally produced state X, see for example [4], [7] - [11]. This is particularly true if 
the outgoing proton momenta can be measured by forward proton taggers: the measured 
proton distributions give spin-parity information about the state X and probe the models 
of soft diffraction used to calculate the rapidity gap survival factors, S 2 [7, 12]. 

Recently there has been much activity in the study of central diffractive processes both 
theoretically (for recent references see, for example, [13]-[16]) and experimentally at the 
Tevatron [17]-[20], by selecting events with large rapidity gaps separating the centrally 
produced state from the dissociation products of the incoming protons. In particular, cen- 
tral exclusive 77 [17], dijet [18, 20] and Xc [19] production have been successfully observed 
at the Tevatron. As noted in [12, 21], these can serve as standard candle processes with 
which we can check our predictions for new physics CEP at the LHC. Indeed, the observed 
rates of all three CEP processes measured at the Tevatron are in broad agreement with 
theoretical expectations [3, 13, 21, 22, 23], which lends credence to the overall theoretical 
framework and motivates further investigation of new and SM CEP physics at the Teva- 
tron, LHC and, as we shall emphasise here, RHIC. Of special interest is the CEP of heavy 
quarkonia since it could provide important information on the physics of bound states and 
can in particular test the current ideas and methods of QCD, such as effective field theories 
and lattice QCD. 

A new area of experimental studies of CEP with tagged forward protons at c.m.s ener- 
gies up to 500 GeV is now being explored by the STAR Collaboration at RHIC [24]-[26]. A 
capability to trigger on and to measure the outgoing forward protons provides an excellent 
means to extend the physics reach in studying CEP processes in exceptionally clean con- 
ditions. The encouraging preliminary results collected in 2009 during Phase I are already 
available [26] and, hopefully, the large data sample expected from the measurements in 
Phase II [25]-[26] should provide some very interesting exclusive physics results. Motivated 
by this, we discuss in this paper the potential for observing exclusive charmonium (x c j and 
rjc) production at RHIC with tagged forward protons (we note that the proton momentum 
loss £ acceptance covers the mass range Mx ~ £-*/s of these states [25]), paying particular 
attention to the new and interesting information that the forward proton distributions can 
provide. We note that the roman pot (RP) detectors at RHIC are uniquely positioned for 
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observing the CEP of such states with tagged protons, as these measurements will not be 
possible at other hadron-hadron colliders in the foreseeable future. 

This paper is organised as follows. In Section 2 we review the theory of the CEP of a 
massive object X in proton-proton collisions. In Section 3 we present predictions for rates 
and distributions corresponding to XcJ and rj c production, paying particular attention to 
distributions in the azimuthal angle difference between the outgoing protons and in the pro- 
ton transverse momenta. Two- and four-body Xc decays could be experimentally important 
observables, and so in Section 4 we discuss the corresponding non-resonant backgrounds. 
Finally, Section 5 contains a summary of our results and some further observations. 

2. Theory 

The formalism used to calculate the perturbative Xc CEP cross section is explained in 
detail elsewhere [7, 12, 23, 27] and so we will only review the relevant aspects here. The 
amplitude is described by the diagram shown in Fig. 1 where the hard subprocess gg —> X 
is initiated by gluon-gluon fusion and the second t-channel gluon is needed to screen the 
colour flow across the rapidity gap intervals. We can write the 'bare' amplitude in the 
factorised form [13, 28] 

T = tt 2 [ d — ; f g (x 1 ,x' 1 ,Qj,v 2 ;t 1 )f g (x2,x' 2 ,Ql,n 2 ;t 2 ) , (2.1) 

where the / ff 's in (2.1) are the skewed unintegrated gluon densities of the proton: in the 
kinematic region relevant to CEP, they are given in terms of the conventional (integrated) 
densities g(x,Qf). ti is the 4-momentum transfer to proton i and \x is the hard scale of the 
process, taken typically to be of the order of the central particle mass: we use fi = Mx /2 in 
what follows. KA is the colour-averaged, normalised sub-amplitude for the gg —> X process 




Figure 1: The perturbative mechanism for the exclusive process pp — > p + X + p, with the eikonal 
and enhanced survival factors shown symbolically. 
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Here a and b are colour indices, Mx is the central object mass, represents the gg — > X 
vertex and qi ± are the transverse momenta of the incoming gluons, given by 

Qi ± =Q±~ Pi ± , Q2 ± = ~Q± - P2 ± , (2.3) 

where Q± is the momentum transferred round the gluon loop while pi ± are the transverse 
momenta of the outgoing protons. Only one transverse momentum scale is taken into 
account in (2.1) by the prescription 

Qi = min{Q_L, |(Qj_ - Pi ± )|} , 

Q 2 = min{Q ± ,|(Q ± + P2± )|} . (2.4) 
Explicitly, the gg — > Xc Vc vertices are given by [22] 



1 c 



V ° = \/^((<7ix<72 ± )(3M^ - q( ± - qij - 2q 2 ± q 2 2± ) , (2.5) 



6M X 

Vi = -^PirtoAfaMtoJ-faMtoJ 2 )^*? > ( 2 - 6 ) 

V 2 = : ^^(s(q 1± Uq 2± ) a + 2(q 1±q2± ) Pl ,p 2a )e*/ a , (2.7) 

V - = ic v (q 1± x q 2± ) • n , (2.8) 



where the normalisation factors are given by 

1 16vra s (/i 2 ) 



Xc 2^iVc (qm) 2 V 47rM ; 



^p(O), 



Xc 



1 4to 5 (^ 2 ) 1 

c - = m • <2 - 9) 



- (S, 2 „h) /d 2 Pu d 2 P2, ^^ffi >>|2 Sg k (PU ■ P2, ) , (2-10) 



and 05(p)(O) is the 5(P)-wave wavefunction at the origin. The renormahsation scale /i 2 is 
assumed to be of the order of m 2 . The differential cross section at X (= % c , Vc) rapidity 
yx is then 

^L_/c2 \^2„ A 2„_ |r(pi ± ,p 2± ))| 2 

dyx 

where T is given by (2.1) and S 2 ^, is the 'eikonal' survival factor, calculated using a gen- 
eralisation of the 'two-channel eikonal' model for the elastic pp amplitude (see [29] and 
references therein for details). For comparison, we will also consider in Section 3 a sim- 
ple one-channel approach, where the elastic amplitude is described by a single Pomeron 
exchange, and is taken to have the Gaussian form 

A pp (s, k 2 ) = is C*a^(s) exp(-Bk 2 /2) , (2.11) 

where we take C* = 1.3 to account for the possibility of proton excitations (p — > N*) 
in intermediate states [30], and we use = 60 mb and B = 14.7 GeV -2 at ^/s = 500 
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GeV [31]. 

Besides the effect of eikonal screening S^, there is some suppression caused by the 
rescatterings of the intermediate partons (inside the unintegrated gluon distribution f g ). 
This effect is described by the so-called enhanced Reggeon diagrams and usually denoted 
as S 2 nh , see Fig. 1. The value of S 2 nh depends mainly on the transverse momentum of 
the corresponding partons, that is on the argument Q 2 of f g (x, x', Q 2 , /x 2 ) in (2.1), and 
depends only weakly on the p± of the outgoing protons (which formally enters only at 
NLO). While S 2 nh was previously calculated using the formalism of [32], we now use a 
newer version of this model [33] which includes the continuous dependence on Q 2 and not 
only three 'Pomeron components' with different 'mean' Q{. We therefore include the 5 en h 
factor inside the integral (2.1), with (S 2 nh ) being its average value integrated over Q±. 

The expected Xc cross section and final-state particle distributions (in particular of 
the outgoing protons) are therefore determined by a non-trivial convolution of the hard 
amplitude T and the soft survival factor S 2 ik . This is modelled in the SuperCHIC Monte 
Carlo [34], which can now explicitly perform the integral (2.10) on an event-by-event basis. 
This allows for an exact generation of the predicted distributions of the final-state central 
particles and outgoing protons, as well as a precise evaluation of the expected cross sections 
after experimental cuts have been imposed, which will depend on both the specific hard 
process gg — > X as well as the effect of secondary rescatterings. 

We end this section with a brief review of the most important uncertainties that are 
present in our calculation, for more details see [12, 13, 22]. Firstly, the cross section depends 
on the conventional gluon densities g(x,Q 2 ) to the fourth power in a region of low x and 
Q 2 , where they are poorly determined. Secondly, we have the uncertainty in the values 
and energy dependence of the non-perturbative survival factors S 2 ik , S 2 nh , which we discuss 
in part in Section 3. Thirdly, the infrared stability of the Q± integral (2.1) depends on 
the presence of the hard mass scale Mx, and it is not completely clear that the x mass 
is large enough to guarantee this. A fourth related uncertainty comes from the possibility 
of a sizeable 'non-perturbative' contribution to the \cj and r\ c rate, where the exchanged 
Pomerons couple directly to the c/c quarks of the charmonium state. In total, if we take 
for each of the four main sources an uncertainty of order ~ * (2 — 3), than adding the errors 
in quadrature we roughly estimate the uncertainty in the total cross section prediction to 
be of order ~ *5. 

However, we would argue that, as our calculation gives a predicted value for the Teva- 
tron cross section that it is in good agreement with the data, it should also give fairly 
reliable estimates for the heavy quarkonium CEP cross sections at RHIC (and LHC) en- 
ergies. Moreover, observables such as the ratios of the predicted perturbative Xcj and r\ c 
cross sections depend weakly on the PDF set used and so carry smaller overall uncertainties, 
although in particular the possibility of sizeable spin-parity dependent nonperturbative con- 
tributions means it is difficult to quantify this statement. The perturbative contributions 
will also receive spin and parity dependent higher-order corrections, which are unknown, 
see [12]. In particular, Xci CEP, which we recall vanishes for on-mass-shell fusing gluons 
(that is at LO) can only occur at NLO, where we account for the fusion of off-shell glu- 
ons. Within the CEP formalism, which makes use of k± factorisation via the explicitly 
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qi ± dependent skewed PDFs in (2.1), we can in part include this initial gluon off-shellness 
at LO to give a non-vanishing, but heavily suppressed, Xd contribution. However, due to 
this suppression, the Xd CEP cross section is highly sensitive to how this off-shellness is 
included in the gg — > Xd matrix element and it may in general also be highly sensitive 
to NLO effects, which can also generate off-shell incoming gluons. To be safe, as in [12] 
we only keep the leading term in the gluon off-shellness qf ± /M 2 to calculate the Xcj and 
rjc CEP rates, but keeping all terms can suppress our predicted cr(xci)/p(Xco) ratio by a 
factor of ~ 3 (the relative Xc2 to XcO ra te is roughly unchanged) 1 . In the case of the lower 
mass 7] c production, this can suppress the (j{j] c ) / 'cr(Xco) ratio by roughly a factor of 2. 

Furthermore, the predictions for outgoing proton distributions are also less affected 
by the uncertainties outlined above, although possible non-perturbative contributions will 
in general affect the predicted distributions, which we calculate only within the pertur- 
bative CEP model. However, we note that in the case of Xc(o,i) an d Vc production the 
proton and (j) dependence of the perturbative and non-perturbative unscreened ampli- 
tudes are approximately equal (up to corrections of order (p±)/(Q'±)), as these follow from 
general symmetry principles [7, 22]. For these states this is therefore not expected to be 
an important source of uncertainty, although in the Xc2 case a large difference between 
the non-perturbative and perturbative distributions cannot be ruled out. In [22] it was in 
particular found that the perturbative and non-perturbative contributions to the total XcO 
cross section are approximately equal c l<r non P ert - « d<7 pert ' (see [13] for details of the model 
used). Previous calculations [36, 37] suggest that the non-perturbative contributions of the 
three Xc states exhibit a similar hierarchy to the perturbative case and so for simplicity 
we may assume, as in [12], the same relative non-perturbative contribution to the Xc(i,2) 
cross section as in the XcO case, and for consistency we make the same assumption for rj c 
production. 

Finally, we note that while the gg — > Xc2 vertex vanishes for J z = fusing gluons 
in the non-relativistic quarkonium approximation, this will in general receive relativistic 
corrections which will allow a J z = component in the production amplitude [22, 23]. 
Although these corrections are expected to be numerically small (see for instance [38, 39]) 
we recall that, due to the J z = selection rule, CEP of non-relativistic Xc2 (in the \J Z \ = 2 
state) is strongly suppressed ((&(x C 2)/c r (Xco) < 1/40)), and therefore these corrections may 
have an effect on the Xc2 rate and distributions. 

3. Cross sections and proton distributions 

In Table 1 we show the predicted differential cross section for Xcj and rj c CEP at yfs = 500 
GeV, calculated following the formalism outlined in Section 2 and [12]. We note that 
the RHIC Xd rate is not expected to be significantly lower than the Tevatron prediction 
of 35 nb (which is in good agreement with the CDF measurement [19] once the Xc(i,2) 
contributions have been included). This is due to the 'eikonal' and 'enhanced' survival 
factors (^eik'^enh)' which both increase with decreasing y^s, and therefore compensate the 
decrease in cross section coming from the smaller gluon density xg(x, Q 2 ) probed at RHIC 

lr The issue of gluon off-shellness in the case of XcO CEP was first addressed in [35]. 
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XcO 


Xcl 


Xc2 


Vc 




27 


0.55 


0.34 


0.24 


^(pp^PP(J/i> + 7)) 


0.31 


0.19 


0.067 





Table 1: Differential cross section (in nb) at y — for central exclusive \c.j and r\ c production at 
\fs = 500 GeV. Also shown is the cross section via the %c — > J/^pJ decay channel. 







XcO 


Xcl 


Xc2 


Vc 


Fraction (%) 


p i± < 0.4 GeV 


47 


17 


19 


13 


Pi ± > 0.5 GeV 


8.6 


9.9 


16 


20 


0.5 < p 1± < 0.7 GeV, p 2± > 0.8 GeV 


1.2 


1.0 


2.2 


1.9 



Table 2: MC cross section fractions for Xcj and ?7 C CEP at y — and \/s = 500 GeV, for different 
cuts on the outgoing proton momenta. S^ ik is calculated using parameter set 1 of the two-channel 
cikonal model of [31], which accounts for the first N* resonance excitation in low mass (p — > N*) 
proton dissociation. 

energies. Indeed, the predicted rj c rate at RHIC is in fact higher than the Tevatron value 
of 200 pb, although the cross section difference is well within the theoretical uncertainties. 

As was shown in [12, 22], the Xc(i,2) rates are expected to be heavily suppressed relative 
to the XcO rate, due to the near-exact J z = selection rule that operates for CEP [7, 23], 
although this suppression may be compensated by the larger Xc(i,2) ~~ > J li'l branching ra- 
tios if Xc CEP is observed via this decay channel. While the lower mass rj c state of course 
does not decay to J/ip'f, it may be observable via 3- or 4-body hadronic decays, which we 
discuss below. In Table 1 we therefore also show the predicted differential cross section for 
Xc CEP at y/s = 500 GeV via the Xc — > J Wt decay channel, which (with photon recon- 
struction using a converter) would represent a promising way to observe the CEP of these 
higher spin states. We show in Table 2 the Xcj and ij c cross section fractions after realistic 
experimental cuts have been imposed on the outgoing protons 2 , to take into account the 
p± acceptance of the roman pot (RP) detectors. In particular, rows 1 (pi ± < 0.4 GeV) 
and 2 (pi ± > 0.5 GeV) correspond to the RPs for Phases I (currently in place) and II 
(to be installed for ~2013) of the STAR pp2pp physics programme, respectively [26]. As 
expected, we can see that by selecting events with small p±, the ratio of the expected XcO 
to Xc(i,2) an d Vc yields increases by a factor of 2-3: as the outgoing proton p± decreases, 
the J z = selection rule becomes more exact and the higher spin states are more heavily 
suppressed. For higher proton p± the effect of absorptive corrections becomes more impor- 
tant: although the 'bare' unscreened Xc(i,2) amplitudes favour higher p± values, this effect 
is less significant once screening effects are included. In particular, the XcO state is not 
significantly suppressed when we select events with high p± , although some enhancement 
of the Xc2 rate is found. On the other hand, in the case of the i] c , for which we recall 
that the unscreened amplitude also favours higher p± values, there is over a factor of 2 

2 Note that the machine constraints may also in general affect the 2iv coverage for the outgoing protons: 
in particular for RHIC Phase II these could give a further factor of ~ 1/(0. 75) 2 decrease in the observed Xc 
rate [26]. 
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Figure 2: Normalised distributions (in arbitrary units) of the difference in azimuthal angle between 
the outgoing protons for the CEP of x c (o.i,2) and r] c states at y = and = 500 GeV, and for a 
range of cuts on the proton p± . S^ ik is calculated using set 1 of the two-channel eikonal model of [31] , 
which accounts for the first N* resonance excitation in low mass (p —¥ N*) proton dissociation. 



enhancement relative to the XcO- 

In Refs. [7, 12], it was shown that the distributions in p± and <f> of the outgoing protons 
depend sensitively on the spin and parity of the centrally produced object. This is illus- 
trated in Fig. 2, where we plot the expected da/d(ft distributions for Xcj and r] c production, 
calculated using the SuperCHIC MC 3 . By applying different cuts to the outgoing proton 
p±, we can also in principle probe the underlying theory in a more detailed way, and to 
illustrate this we plot the <f> distributions for three different sets of p± cuts: P(i,2) ± < 0-4 
GeV, P(i,2) ± > 0.6 GeV and 0.5 < pi ± < 0.7 GeV, p2 ± > 0.8 GeV. For low p± the screening 

3 Note that the unscreened Xcl an d Xc2 distributions, shown in [12], are different from those presented 
in the published version of [15]. However we have been informed by the authors that their more recent 
calculation confirms our results [40]. 
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XcO 


Xcl 


Xc2 




Fit 1 


0.092 


0.23 


0.15 


0.26 


w> Fit 2 


0.062 


0.18 


0.11 


0.21 


(S^ik)' Single channel 


0.070 


0.20 


0.13 


0.24 



Table 3: 'Eikonal' survival factor, (S 2 ik ), averaged over the outgoing proton pj_, calculated using 
the two-channel eikonal model of [31], with two different choices of model parameters. Parameter 
set 1 accounts for the first N* resonance excitation in low mass (p — > N*) proton dissociation, 
while set 2 includes excitations up to a larger M 2 ~ 6 GeV 2 . Also shown is the result of using the 
simplified single channel eikonal approach, see (2.11). 
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Figure 3: Normalised distributions (in arbitrary units) of the difference in azimuthal angle between 
the outgoing protons for \c0 CEP at y = 0, with the survival factor S 2 ik (bt) calculated using the 
two-channel eikonal model of [31], with two different choices of model parameters. Parameter set 1 
accounts for the first N* resonance excitation in low mass (p —} N*) proton dissociation, while set 
2 includes excitations up to a larger M 2 ~ 6 GeV 2 . Also shown is the result of using the simplified 
single channel eikonal approach, see (2.11). Note that the 'two-channel, set 1' XcO distributions are 
the same as those plotted in Fig. 2. 



corrections do not change this 'bare' behaviour too much, however in the case of a rela- 
tively large p± (green and blue lines in Fig. 2) the role of absorptive effects becomes quite 
visible: starting from = the absorptive correction increases with <f> producing a dip in 
the region of 4> ~ 7r /2 for the cases of the XcO and Xc2 and about <p ~ 2.3 for the Xci- We 
note that these characteristic 'diffr active dip' structures have the same physical origin as 
the proton azimuthal distribution patterns first discussed in [29]. For the r\ c these effects 
are less significant, although some non- negligible dip structure around <f> = 2 — 2.5 can be 
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0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 



pi (GeV 2 ) 

Figure 4: Normalised distributions (in arbitrary units) of the p\ of the outgoing proton, integrated 
over the other proton p±, for \c.j and rj c CEP at y = 0, and for \c0 CEP with the survival factor 
S 2 ik (b t ) calculated using the two-channel eikonal model of [31], with two different choices of model 
parameters. Parameter set 1 accounts for the first N* resonance excitation in low mass (p — > N*) 
proton dissociation, while set 2 includes excitations up to a larger M 2 ~ 6 GeV 2 . Also shown is the 
result of using the simplified single channel eikonal approach, see Eq. (2.11). The lower Xcj plot 
uses fit 1. 



seen. 

As described in [22], a further way to extract spin information about the centrally pro- 
duced x state is by measuring the angular distributions of its decays products, in particular 
here the final state pair from J/tf) decay. These spin-dependent angular distribu- 

tions, which are generated by the SuperCHIC MC, would also represent an interesting 
observable, providing complementary information to the tagged proton distributions. 

As well as depending on the spin-parity of the centrally produced particle, the proton 
distributions will also be affected non-trivially by soft-survival effects [29], in particular 
the p_i_ dependent survival factor S^jj^pij., P2±)> through Eq. (2.10). We show in Fig. 3 
the predicted <f> distribution for XcO CEP using the 'two-channel' eikonal formalism of [31] 
for a range of cuts on the proton p±. To give some indication of the uncertainty involved 
in calculating the survival factors, we show the distributions using two different choices 
of model parameters as well as using the simplified single channel model described in 
Section 2. In particular, to calculate S^ ik the elastic and cr tot pp data were fitted and the 
available data on low mass (p —> N*) proton dissociation were used, as described in [31, 41]. 
Parameter set 1 and the one-channel enhanced (by C = 1.3) eikonal model account for the 
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fixed target (mainly old FNAL) data and correspond to the first N* resonance excitation. 
Set 2 accounts for CERN-ISR data and excitations up to a larger M 2 ~ 6 GeV 2 (which is 
correlated with S 2 nh based on the KMR model [42]). 

In Table 3 we show the averaged suppression factors (S 2 ik ), for the different param- 
eter sets and the simplified single-channel model. Clearly there is some (of order < 50%) 
non-negligible model dependence in the overall soft suppression factor, which varies de- 
pending on the spin-parity of the centrally produced particle. This variation, which results 
from our incomplete knowledge of the soft amplitude used to calculate the survival factors, 
represents an important source of uncertainty in any CEP cross-section predictions. On 
the other hand, the ratios of the S 2 ik values between the different Xcj and r\ c states are less 
model dependent, varying by < 15%. Moreover, the distributions are not too dependent 
on the choice of parameter set, although some difference, in particular for larger proton p± 
values (where we recall the amplitude is more sensitive to soft survival effects), is apparent. 
This is equally true for the case of the single-channel eikonal model, in particular when 
experimental cuts are imposed. We note that the Xd and Xc2 4> distributions are also 
similarly dependent on the choice of parameter set for the survival factors. 

Finally, we show in Fig. 4 the proton p\ distribution (integrated over the other proton 
p±) for XcO CEP using the same two sets of parameters and the simplified single-channel 
model to calculate the survival factors. When no cuts are imposed on the proton pj_ we 
can see that there is in principle an observable difference between the p\ distributions 
resulting from the different parameter sets, but this is no longer the case when the most 
relevant cut pi ± > 0.5 GeV is imposed. As with the 4> distributions, we can see that the 
p\ distribution (shown for the case pi ± > 0.5 GeV) also depends on the spin-parity of the 
centrally produced state. 

XcO CEP, observed via two-body (e.g. XcO — > tt + tt~, K + K~, pp) or four-body (e.g. 
XcO — > 2(-7r + 7r~), 7r + 7r~ K + K~) decays also represents an interesting observable at RHIC, 
provided the direct QCD backgrounds are sufficiently under control — we will discuss this 
below. Experimentally, the exceptionally good resonance mass reconstruction (of order of 
a few MeV) provided by the excellent charge particle identification and high resolution 
tracking in the STAR TPC will greatly aid in increasing the S/B ratio. On the other 
hand, the higher spin Xc(i,2) states are expected to give negligible contributions via these 
decay channels. To give an indication of the expected rates, in Table 4 we summarise the 
branching ratios, taken from [43], for some of the two- and four-body XcO decays that may 
be experimentally relevant at RHIC. Note that around 6% of all XcO decays are to h + h~ or 
2{h + h~) with h = ir,K. The corresponding mesonic two-body decays (tttt, KK, cfxp, ww, r]t] 
etc.) of the Xc2 are all at the 10 -3 level, while the 4-body decay modes (e.g. 7r + 7r~ K + K~ , 
2(7r + 7r~)) are at the 1% level [43], and may also be useful. 4 In the case of r\ c production, 

4 We note that currently the STAR mass acceptance for 7r + 7r~ and K + K~ states does not extend 
beyond ~ 3.2 GeV due to particle ID, i.e. below the \ c threshold [26]. However, even without particle ID, 
observation of XcO — > tt + tt~ , K + K~ should be possible: in particular, for each event we can include the 
possibility in the analysis that the two charged tracks are either it + it~ or K + K~ . This will increase the 
non-resonant background by about a factor of 2, but given the excellent mass resolution in the STAR TPC, 
a clear \c peak should still be observable. Nevertheless, any future extension of the STAR mass reach (with 
particle ID) would certainly improve the experimental situation. Note also that the mass coverage of the 
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Mode 


Branching ratio 


7T + 7r~ 


(0.56 ±0.03) x 10~ 2 


K+K- 


(0.610 ±0.035) x 10~ 2 


K S K S 


(0.316 ±0.018) x 10~ 2 


H> 


(9.2 ± 1.9) x 10" 4 


fo (980) f (980) 


(6.8 ±2.2) x 10~ 4 


pp 


(2.28 ±0.13) x 10~ 4 


AA 


(3.3 ± 0.4) x 10" 4 


2(7r+vr-) 


(2.27 ±0.19) x 10~ 2 




(1.80 ±0.15) x 10" 2 



Table 4: Branching ratios for Xco two- and four-body decays, taken from Ref. [43]. 



the three-body (e.g. KKir, with branching 7%) and four-body (e.g. direct 2(7r + 7r~) decay, 
with branching 1.2%, or via the wide resonance states p°p° — > 2(7r + 7r~), with branching 
2%) decay modes are the most realistic. 

4. Non-resonant QCD background 

In Section 3 we noted the possibility of observing Xco CEP via two-body decay channels, 
with the Xc ttt decay being a promising example. We recall that these decay channels, 
especially 7T7t, K + K~ and pp, are ideally suited for spin-parity analysis of the Xc states, in 
particular the fact that the Xc(i,2) t wo body branching ratios are in general of the same size 
or smaller (or even absent for the Xd) than the Xco ensures that the J z = selection rule 
is active, see [13, 22] for more details. However, in this case we may in principle expect a 
sizeable background resulting from direct production, pp — ^ p + hh + p with h = ir, K,p, 
and so care must be taken to estimate the expected non-resonant contribution. This 
can be modelled using a 'non-perturbative' framework, mediated by Pomeron-Pomeron 
fusion with an intermediate off-shell h = ir, K,p exchanged between the final-state particle 
pair, see Fig. 5(a) 5 , or via the perturbative 2-gluon exchange mechanism with the gg — > 
tt + tt~,K + K~ coupling modelled using the formalism of [47, 48] 6 , see Fig. 5(b). In the 
kinematic (M nn ~ M x , p±(ir) ~ M x /2) region relevant to Xc production we may expect 
both mechanisms to contribute to the overall cross section. Besides the 'bare' amplitudes 
shown in Fig. 5 we also have to account for the absorptive corrections, that is the eikonal, 
5 e ik, and 'enhanced', S^h, survival factors introduced in Section 2. 

Considering first the perturbative mechanism for 7r + 7r~ production, using purely di- 
mensional arguments we can see that the amplitude to form an exclusive pion pair with 
large transverse momentum k± will be proportional to the ratio {f-w/k^) 2 (the pion form 
factor f n enters through the normalisation of the pion wavefunction (f>(x)), that is the 

pp decay mode is expected to be broader (up to around 4 GeV) due to better particle ID. 

5 7T7r CEP mediated by double Pomeron exchange has been a subject of theoretical studies in the last 40 
years or so (see [44, 45] for early references and [46] for a more recent one). There have also been a variety 
of experimental results, in particular from the CERN ISR (for recent reviews see [5, 6]). 

6 We note that gg — ¥ pp can also in principle be modeled in this way, see [49, 50]. 
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(a) 



(b) 



Figure 5: (a) Non-perturbative tttt CEP production mechanism - similar diagrams exist for 
K + K~ and pp production, (b) A typical diagram corresponding to the perturbative gg — > tttt 
hard subprocess. 

cross section of the gg pp — > tt + tt~ hard subprocess contains the numerically small fac- 
tor (iV /£;_!_ ) 4 . Moreover, as we shall discuss in detail in a future publication [51], the LO 
amplitude for tttt exclusive production in a J z = state vanishes in the same way as the 
77 —7- 7T°7r amplitude [47, 48] . Thus the process will occur mainly through the admixture 
of the | J z \ =2 state of the incoming active gluons gg pp , which we recall is strongly sup- 
pressed due to the J z = selection rule which operates for forward outgoing protons [23]. 
We note that this result also extends to K + K~ production by replacing f w by fx and 
using a slightly different (asymmetric) wavefunction to account for the non-zero strange 
quark mass [48]. The net effect is that the K + K~ cross section is predicted to be approx- 
imately equal to that for tt + tt~ (this has already been observed by BELLE for the case of 
77 7T+TT-, K + K~ [52], see also [53]). 

In the case of the non-perturbative contribution, large pion p± values are suppressed 
by the form factor, OI the intermediate off-shell pion, although in this case it is not 

completely clear what form to take for F n (t); as a result there is a large uncertainty in the 
predicted cross section for this non-perturbative mechanism. In particular, we may either 
take the 'soft' exponential exp(fti), or the 'hard' power-like form ~ 1/t, both of which are 
used in the literature [45, 46, 54]. As \t\ (that is, the pion or kaon p±) is increased, these 
two choices of form factor give vastly different cross sections: while the exponential form 
gives a negligible cross section, the power-like form gives a comparable cross section to Xc 
CEP. When available, CDF data [6, 55] on central exclusive tt + tt^ production in the Xc 
mass range, collected in the same rapidity range (\rj\ ~ 1), would therefore help to further 
constrain the non-resonance background and may provide useful information about the off- 
shell pion form factor. However, we note that even taking the power-like pion form factor, 
which gives the much larger non-resonant cross section estimate, we expect to obtain a 
background contribution that will be under good control experimentally. In particular, it 
corresponds to a cross section of about 10-20 pb under the Xc resonance mass peak — an 
order of magnitude less than the expected XcO — > 7r + 7r~ CEP cross section. This is equally 
true for the K + K~ channel. 

For the pp decay channel, we again have to consider the production mechanism shown 
in Fig. 5 (a), but with a proton, rather than pion, line. The Pomeron-proton coupling 
is about a factor 3/2 larger than in the pion case, but the proton form factor is steeper 
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(the 'hard' power-like form is given by F p ~ 1/i 2 )- Therefore, for a reasonably large 
proton k± > 1 GeV the non-perturbative contribution is expected to be smaller than in 
the 7T7r-channel. 

In the case of 4-body decays, the situation is more complicated. An important part of 
these signal decays come from the Xc decay to two wide resonances, such as Xc — > P°P° — > 
2(7T + 7r~). An analogous mechanism is possible for the background process, and, taking 
into account the number of appropriate resonances and the final state interactions, it is 
hard to give a definite prediction. However, we do not expect the 4-body background to 
be much larger than for the 2-body channel. 

5. Conclusions 

In this paper we have studied Central Exclusive Production of Xcj and rj c mesons at RHIC, 
for proton-proton collisions. Because the gap survival factors increase with decreasing 
■s/s, the production rates at RHIC are in fact expected to be comparable to those at the 
Tevatron. For the Xcj — > J/^> + 7 decay channel, the strong suppression of the J = 1,2 
states is compensated by the larger branching ratios for these states. Overall, for this decay 
channel we predict a cross section do /dy x ~ 0.6 nb at y x = 0, with the J = 0, 1, 2 states 
roughly in the proportion 5:3: 1. We have also considered rj c CEP, for which we predict 
a cross section of approximately 0.2 nb, which is heavily suppressed relative to the XcO (for 
which the predicted rate is roughly 30 nb), as we expect from the = + selection rule. 
These rates are, however, dependent on the experimental p± cuts on the outgoing protons. 
In particular, selecting events with small p± increases the ratio of the XcO to Xc(i,2)/Vc yields 
by roughly a factor of 2: as the outgoing proton p± decreases, the J z = selection rule 
becomes more exact and the higher spin states are more heavily suppressed. SuperCHIC 
MC [34], which generates all of the processes described in this paper with exact proton 
kinematics, allows for a complete evaluation of the acceptances within the full experimental 
setup for the different cc states. 

The distribution in the azimuthal angle difference between the outgoing protons also 
depends sensitively on the spin and parity of the centrally produced object, and also to 
some extent on the transverse momentum dependence of the survival factors. While we 
focused mainly on the Xcj — > J/ip + j decay channel, we also considered two- and four- 
body decays which are particularly relevant for XcO production. We showed that in the 
two-body case, backgrounds from 'direct' QCD production based on both perturbative and 
non-perturbative models are expected to be under control. 

It is worth mentioning that CEP can help to shed light on the nature of the numer- 
ous recently discovered charmonium-like mesons (X,Y,Z) [43], and, in particular, on the 
X(3872) resonance, whose nature is presently unclear [56]. Found in 2003 by the BELLE 
collaboration [57] and then confirmed by CDF [58] and DO [59] at the Tevatron and by 
BABAR [60], this narrow state still remains a mystery and, consequently, a subject of in- 
tensive theoretical activity (for reviews see for example [61, 62] and for recent papers [63]- 
[67] and references therein). As is the case for many of the new charmonium-like states, 
A (3872) cannot be easily fitted into the conventional charmonium spectroscopy, and so 
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a number of more exotic explanations for its origin have been suggested, such as a DD* 
molecule, tetraquark, hybrid, some sort of threshold effect, etc. — see [61, 62] for more 
details. However, quite recently there has been a renewal in attempts to interpret it as a 
regular l ++ (2 3 i => i) charmonium state [64] - [67]. 7 

Additional experimental information on the properties of the A (3872), especially its 
J p assignment, could provide new insight into the theoretical understanding of this state. 
A detailed study of the CEP process pp — > p + X(3872) + p would be a valuable new 
tool to test the quantum numbers of this state (see [13]), and so could play an important 
role in the resolution of the X(3872) puzzle. In particular, a comparison of the proton 
momentum correlation pattern with the expectations for a 1 ++ charmonium state would 
provide support for or rule out the interpretation of this resonance as a 2 3 Pi member of 
the charmonium family. 

Finally, we note that the planned vertex detector will allow STAR to perform a search 
for the Xc2(2-P) resonance, which has recently been observed at B-factories [69] in the DD 
mode. We do not expect that the CEP rate for this state will be very different, within the 
theoretical uncertainties, from the prediction given here for the % C 2(1P) state. 
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